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ABSTRACT 

A particularly interesting discovery in observations of GRB 121027A is that of a giant X-ray bump 
detected by the Swift/XBT. The X-ray afterglow re-brightens sharply at ~ 10'^ s after the trigger by 
more than two orders of magnitude in less than 200 s. This X-ray bump lasts for more than lO** s. 
It is quite different from typical X-ray flares. In this letter we propose a fall-back accretion model 
to interpret this X-ray bump within the context of the collapse of a massive star for a long duration 
gamma-ray burst. The required fall-back radius '-^ 3.5 x 10^° cm and mass ~ 0.9 — 2.6 Mq imply that 
a significant part of the helium envelope should survive through the mass loss during the last stage 
of the massive progenitor of GRB 121027A. 

Subject headings: accretion, accretion disks-black hole physics-gamma-ray burst: individual (GRB 
121027A)-magnetic fields 
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1. INTRODUCTION 

The most popular models of long-duration gamma- 
ray bursts (GRBs) invoke a collapse of a massive star 
(Woosley 1993; Paczynski 1998; MacFadyen & Woosley 
1999). The accretion of the stellar core by the central 
black hole (BH) fuels the prompt emission. However, the 
central engine activity does not cease after the prompt 
phase, which is especially supported by observations of 
late X-ray flares by the X-Ray Telescope (XRT) onboard 
the Swift satellite. X-ray fiares share a lot of similarities 
with GRB prompt emission and are therefore interpreted 
by the same mechanism as prompt emission, i.e., inter- 
nal shocks (Burrows et al. 2005; Zhang et al. 2006; Fan 
& Wei 2005). The fall-back and accretion of the stel- 
lar envelope may have some observational consequences 
(Kumar et al. 2008a, 2008b; Dai & Liu 2012). In order 
to have a successive GRB jet penetrating the star, the 
progenitor is usually assumed to be a Wolf-Rayet star 
that has an evolved compact Helium envelope. 

For the GRB BH central engine models, there are two 
main energy reservoirs to provide the jet power: the 
accretion energy in the disk which is carried by neu- 
trinos and anti-neutrinos, that annihilate and power a 
bipolar outfiow; and the spin energy of the BH which 
can be tapped by a magnetic field connecting the outer 
world through the Blandford-Znajek (1977, hereafter 
BZ) mechanism. Both models have been extensively in- 
vestigated by many authors (e.g., Popham et al. 1999; 
Lee et al. 2000; Li 2000; Narayan et al. 2001; Di Matteo 
et al. 2002; Kohri & Mineshige 2002; Wang et al. 2002; 
McKinney 2005; Gu et al. 2006; Chen & Beloborodov 
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2007; Janiuk et al. 2007; Lei et al. 2009; Liu et al. 2007). 

GRB 121027A may provide us a good chance to study 
the properties of GRB progenitor as well as the central 
engine models. As presented in the next section, the ris- 
ing behavior of the giant X-ray bump in GRB 121027A 
is quite different with typical X-ray flares, so it should 
have different physical origin. In Section 2, we describe 
the prompt trigger and late XRT observations of GRB 
121027A. In Section 3, we propose the fall-back accretion 
model and apply this model to the giant X-ray bump ob- 
served in GRB 121027A. In Section 4, We briefly sum- 
marize our results and discuss the implication. 

2. GRB 121027A OBSERVATIONS 

GRB 121027A was discovered at Tq = 07 : 32 : 29 
UT on 2012 October 27 by the Burst Alert Telescope 
(BAT) on board Swift (Evans et al. 2012a) and was 
later accurately located by XRT at a position of a = 
04h 14m 23.37s, S = -58°49'46.6" (J2000), with an un- 
certainty of 2. (Beardmore et al 2012). The redshift 
of this burst is provisionally implied to be z ^ 1.77 by 
Gemini South spectroscopic observation, assuming the 
significant absorption feature at 7770A in the optical 
spectrum as the Mgll doublet 2796/2803A (Tanvir et 
al. 2012). It is later confirmed and accurately measured 
to be z = 1.773 by identifying several metal absorption 
lines with the X-shooter spectragraph (Kruehler et al. 
2012). The results of data analysis of the prompt BAT 
observation are as follows (Barthelmy et al. 2012). The 
mask-weighted light curve shows a pulse with fast-rise 
and exponential-decay (FRED) profile plus two small 
peaks on the tafl. The duration of GRB 121027A is 
Tgo = 62.6 ± 4.8 s in the 15 - 350 keV band. The 
time-averaged spectrum is best fit by a simple power 
law with the photon index F = 1.82 ± 0.09. The flu- 
ence in the 15 - 150 keV band is f-y = 2.0 ± 0.1 x 10^^ 
yielding an isotropic gamma-ray energy 



erg cm ^ s ^ 

release E^^i^o = iTrDlf^/{l + z) ^ 1.58 ± 0.08 x 10' 



,52 



erg. Here we adopt the concordance cosmology with 
Qm = 0.27, Qa = 0.73 and ho = 0.71. 

The XRT began observing the burst at Tq + 67.4 s. 
The XRT light curve shows several components (Evans 
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Fig. 1.— BAT (gray) and XRT (black) light curves of GRB 
121027A. BAT flux is calculated at 10 keV. XRT flux is absorption- 
corrected in the 0.3 — 10 keV. The prompt and afterglow emission 
is separated by the dashed vertical line. Also plot is the Japanese 
MAXI/GSC observation (open pentagram), whose flux has been 
extrapolated from the 4 — 10 keV band (3.6 X 10~^ erg cm~^ s~^) 
to the 0.3 - 10 keV band (8.3 X 10"^ erg cm-^ g-i). 

et al. 2012b). FigH] shows the XRT hght curve, which 
is the temporal evolution of the unabsorbed flux in the 
0.3 - 10 keV band (Evans et al. 2007, 2009). Initially it 
decays as a power law with the temporal index ai ^ 1.8 
until an X-ray flare emerges at Tq + 220 s and lasts for 
^ 300 s. The decay slope of this X-ray flare is very steep, 
a2 ~ 6.6. 

The most interesting discovery in GRB 121027A is 
the giant X-ray bump in the subsequent observations. 
The flux of the giant X-ray bump increases sharply at 
~ To -F 10^ s by more than two orders of magnitude in 
less than 200 s. The flux is 1.2 x 10"^-^ erg cm~^ s~^ 
at To -I- 1033 s (Windowed Timing mode), and suddenly 
increases to 1.7 x 10"^ erg cm^^ s^^ at To + 1198 s (Pho- 
ton Counting mode). Such "step-like" re-brightening of 
the X-ray bump in GRB 121027A is quite different with 
those of typical X-ray flares detected by Swift in the 
past eight years. For X-ray flares, the rise and decay 
time scales are compatible, which are close to the peak 
time of the flare (e.g., Liang et al. 2006; Chincarini et 
al. 2010). For the giant X-ray bump of GRB 121027A, 
the decay time scale, as shown in Figure 1, is ^ 10^ 
s, much longer than the rising timescale. There is no 
XRT observation between To + 1.2 ks and Tq + 5.3 ks; 
however, the MAXI/Gas Slit Camera(GSC) detects the 
X-ray counterpart of GRB 121027A with the flux 150 
mCrab in the 4-10 keV band at To + 2424 s (Serino 
et al. 2012). Although the decay of the bump is not 
smooth and may have flaring features, its envelope can 
be roughly divided into two stages. Initially the decay 
slope is 1.6 until To + 1.2 x lO'' s. After that, the de- 
cay slope is ~ 3.8. The spectrum gets softening during 
the decay phase. The time-averaged photon index F is 
1.66 ± 0.20 between To -I- 5.3 ks and Tq + 11.9 ks, while it 
changes to a much softened value of 2.47 ± 0.13 between 
To + 16.3 ks and To + 30 ks. The late XRT emission is 
dominated by another component after To-|-50 ks, which 
can be attributed to the forward shock emission. More 
detailed data analysis of spectra and lightcurves of BAT 
and XRT emission of GRB 121027A can be found in Peng 



et al. (2013) and Levan et al. (2013). 

3. FALL-BACK ACCRETION MODEL 
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Fig. 2. — Illustration of our model. The fall-back and accre- 
tion of the stellar envelope produce the X-ray bump seen in GRB 
121027A. The jet is powered by the Blandford-Znajek mechanism, 
which extracts the rotational energy of Kerr BH through large-scale 
magnetic field (oblique and spiral lines). 



We suggest that the X-ray bump seen in GRB 121027A 
is the result of the fall-back accretion, as shown in Figl^J 
The fall-back accretion is expected to start at the fall- 
back time ^0 = tfh (in the following, time is defined in 
the cosmologically local frame). The fall-back time is 
the time it takes a parcel of gas of the progenitor star at 
radius rfb to fall to the center, and it is approximately 
equal to the free-faU time, ti^ ~ {n'^rf^/SGM,)^/'^ , where 
M, is the BH mass. 

Following MacFadyen et al. (2001), Zhang et al. 
(2008), and Dai & Liu (2012), the fall-back accretion rate 
initially increases with time as Mcariy oc t^^"^ until it 
reaches a peak value at t-p. The late-time fall-back ac- 
cretion behavior follows Miatc oc t~^/^, as suggested by 
Chevalier (1989). Therefore, we assume that the fall- 
back accretion rate evolves as a smooth-broken-power- 
law function of time|3 



M = Afn 



tp — to 



t-tn 



to 



1 



(1) 

where a,. — 1/2, ad — —5/3 and s describes the sharp- 
ness of the peak. The dimensionless accretion rate is 
defined as m = M/{Mq/s). 

As suggested by Lei et al. (2013; see also Lei et al. 
in preparation), the jet may be dominated by the BZ 
power especially at late time (Fan et al. 2005; Zhang 

The actual peak time and peak accretion rate in Eq. (1) are 
slightly different from tp and Mp. 
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3 



& Yan 2011). For this reason, we connect the observed 
X-ray luminosity to the BZ power through 

vEb = fhLx,iso (2) 

where 77 is the efficiency of converting BZ power to X-ray 
radiation and /b is the beaming factor of the jet. 

The BZ jet power from a BH with mass M, (or dimen- 
sionless mass m, = M,/Mq) and angular momentum J, 
is (Lee et al. 2000; Li 2000; Wang et al. 2002; McKinney 
2005; Lei & Zhang 2011; Lei et al. 2013) 



Eb = 1.7 X lO^^a^m^B; 15^(0.) erg s" 



(3) 



where 5.^5 = BJW^G and 

Fia,) = [(1 + q^)/q^][{q + 1/q) arctanq - 1]. (4) 

Here a, = J,c/ (GM^) is the BH spin parameter, and q = 
a,/{l + al). For < a, < 1, 2/3 < F{a,) < 7r-2. 

It is obvious that the BZ power depends on M,, B,, and 
a,. A strong magnetic field of ~ lO^^G is required to 
produce the high luminosity of a GRB. 

As the magnetic field on the BH is supported by the 
surrounding disk, there are some relations between B, 
and M. In a hyper-accreting flow in a GRB, it is possible 
that a magnetic flux is accumulated near the BH horizon. 
Considering the balance between the magnetic pressure 
on the horizon and the ram pressure of the innermost 
part of the accretion flow (e.g. Moderski et al. 1997), 
one can estimate the magnetic field strength threading 
the BH horizon 



Svr 



pc 



Mc 



(5) 



where r, = (1 + -^/l — a^)?'g is the radius of the BH 
horizon, and — GM,/(?. It can be rewritten as 



B. ~ 7.4 X 10iSmi/2m7i(l + alY^G. (6) 

Inserting it to Equation ([3]), we obtain the magnetic 
power as a function of mass accretion rate and BH spin, 
i.e. 

Eb = 9.3 X lQ^^aimX{a,) erg s~\ (7) 

and 



X{a,)^F{a,)/{l + ^l-alf 



(8) 



It is found that X(0) = 1/6, and X{1) = tt - 2. In gen- 
eral, a faster BH is more favorable for GRB production, 
as w revealed by recent GRMHD numerical simulations 
(Nagataki 2009). 

The BH evolves with time during a GRB, since it would 
be spun up by accretion and spun-down by the BZ mech- 
anism. The evolution equations of a Kerr BH in the BZ 
model can be written as 

^Mc^E^,-Eb, (9) 



dt 

dJ, 

IT 



ML, 



Tb 



(10) 



Since a, = J,c/{GM^), by incorporating Eqs. (|9]) and 
pI7|) . the evolution of the BH spin can be expressed by 



da, , ■ 

-dt = 



Tb)c/(GM2) 



2a.{Mc^E^, - Eb)/{M,^), (11) 

where £'i„s and L,ns are the specific energy and angular 
momentum corresponding to the inner most radius rms 
of the disk, respectively, which are defined as (Novikov 
& Thorne 1973) 



Ems — 



4\/-Rms - 3a. 

%/3-Rms 



GM,2{i^/R^,-2a,) 



V3/R; 



(12) 



(13) 



where R-aia = f^s/fg is the radius of the marginally stable 
orbit in terms of Vg. The expression for R^g is (Bardeen 
et al. 1972), 



R„ 



3 + Z2-[(3-Zi)(3 + Zi + 2Z2)]'/', (14) 

for < a. < 1, where Zi = 1 + (1 - a^)^/^[(l + a,y/^ + 
(l-a.)i/3],Z2^(3a2 + Z2)i/2. 

In Eq. (jlOp , Tb = Eb / is the total magnetic torque 
apphed on the BH, i.e. 



Tb = 3.36 X lO^^^^^'m^B^ 15^(0.) g cm^ s"^ 

(15) 

Here — 0.5il, is usually taken to maximize the BZ 
power, and 



a. 



GM. 2(1 + yr 



(16) 



is the angular velocity of the BH horizon. 

One can calculate the evolution of the BH spin by 
inserting the above expressions into Eq.(ll), and then 
study the time-dependent X-ray luminosity by substitut- 
ing the evolutions of the BH sin and fallback accretion 
rate into Eq.Q of the BZ power. 

For GRB 121027A, the X-ray bump appears at 
To + 1000 s after the GRB trigger, which, divided by 
1 + z, corresponds to tfh ^ to ^ 360 s. This suggests the 
minimum radius around which matter starts to fall back 
is 

Tfb ~ 3.5 X 10^°(Af./3MQ)i/3(tfb/360 sf/^cm. (17) 

The peak flux of the bump in the 1 — 10 keV band is 
Tx.pcak > 8.3 X 10~^ erg cm^^ s~^, corresponding to a 
peak luminosity TxTso > 1-8 x lO'^" erg s"^ From Eq. 
^ , we can estimate the peak accretion rate 

Mp ~ 1.1 X lO~^LxMo,5oa-^X-\a.)7]Zlh.-2MQ s-\ 

(18) 

where ?7_2 = ??/10 ^ and /b -2 = /b/10 ^. 

The total rising time of the bump is about 1800 s (Fig. 
3), we thus have ip - to - 1800/(1 z) s - 650 s. By 
Eq.([T|), the total fallback/ accreted mass should be 



fb 



Mdt - 2Mp(ip -to)/ 3 



to 



4.6 X 10^^Lx.iso..5oa,^X-\a,)7jZlfb,-2Mr, 



(19) 
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Fig. 3.— Fitting (solid line) to the X-ray bump of GRB 121027A 
(a, = 0.9). The inset panel shows the evolution of the dimension- 
less mass accretion rate. 

From Eq.(|6l), the maximum magnetic field strength 
around BH is 

(20) 

The exact values of the above parameters depend 
strongly on the BH spin a, at time tp. 

To constrain the above parameter values in GRB 
121027A, we carried out numerical calculation of Eqs.(IT]) 
- ((TT|) . We obtain the time evolution of the BZ power, 
and compare it with the observations of the X-ray bump 
in GRB 121027A. In our calculation, we assume rj = 
10^^, /b = 10^^ (the jet half-opening angle is con- 
strained to be 0j > 0.2 radian by late XRT observation; 
Peng et al. 2013). The BH is initially set up with a mass 
m, = 3 and a spin a, = 0.9. The calculation starts at 
^0 — 1150/(1 -I- z) s. FiglHshows our model fit to the X- 
ray bump in GRB 121027A. The parameters for the fit- 
ting are nip = 6.1 x 10^**, s = 1.9, and tp = 2950/ (I + 2) s. 
The total fall-back mass is Mfb = 0.9Mq. The total fall- 
back mass is one order of magnitude higher than that es- 
timated with Eq. (19), which underestimates the actual 
duration of smoothed peak accretion (see Fig. 3). The 
initial spin of the black hole for the X-ray bump phase 
has large uncertainty. We also considered low spin cases. 
For a, = 0, the model parameters are rhp — 7.0 x 10"'^, 
s = 0.35, tp = 1350/(1 + z) s, and Mfb = 2.6Mq. For 
a, = 0.5, the model parameters are nip = 1.7 x 10"'^, 
s = 0.70, tp = 2050/(1 + z) s, and A/fb = I.SMq. 
Note that for the decay phase of the bump, we focus on 
the temporal evolution of the envelope of XRT emission. 
There may be fragmentation during the fall-back phase 
(King et al. 2005), which can account for the variations 
of XRT flux during the decay phase. 

There is a break in the lightcurve at time tb ^ 1-2 x 
lO^s, after which the photon index F changes from 1.66± 
0.20 to 2.47 ± 0.13. If the central engine ceases at this 
time, one can expect a transition in the lightcurve from 
the fallback phase to the tail emission phase. In this 
scenario, the temporal decay index after ^b is a = 1 + 



F = 3.47 ± 0.13 due to the curvature effect, which is 
consistent with the observed one (3.8). It requires that 
the mass fallback should stop at ib- Therefore, remnant 
emission comes from high latitude and the observed XRT 
flux fades as t~^-^ for t > t\j. 

4. CONCLUSIONS AND DISCUSSION 

The X-ray afterglow light curve of GRB 121027A is 
unusual. The "step-like" re-brightening at about 1000 s 
since the burst with a duration longer than 10^ s, which 
we refer to as the giant X-ray bump in this Letter, is quite 
different with typical X-ray flares observed by Swift. We 
propose a fall-back accretion model to interpret this X- 
ray bump within the context of the collapse of a mas- 
sive star for long duration GRBs. The fallback radius 
of rfb ~ 3.5 X 10^" cm and mass Mfb ~ 0.9 - 2.6 M© 
for this burst require the helium envelope of the pro- 
genitor should be partly survived before the ending of 
the massive star. One may ask why this burst shows 
the fall-back signature, while most other long GRBs do 
not. One should always have fall-back. The reason may 
be that in the collapsar models, the bounding shock re- 
sponsible for the associated supernova transfers kinetic 
energy to the envelope materials. The more energetic the 
supernova shock, the less envelope material falls back 
into the center. The potential energy of the fall-back 
material at rfb is negative, with an absolute value is 
GM.Mfb/rfb ~ 2x 10'*9(Mfb/Mo) erg for GRB 121027A, 
assuming M, ^ 3Mo. If the kinetic energy delivered 
from the supernova shock is less than the potential en- 
ergy, then these material will fall back. In our scenario, 
GRB 121027A might be accompanied with a low energy 
supernova, or even a failed supernova. 

In this work, we do not include the magnetic coupling 
effect (see Eqs. l9l and fTO|) between the BH and the disk 
through closed magnetic field lines (Li & Paczynski 2000; 
Wang et al. 2002; Lei et al. 2009; Janiuk & Yuan 2010). 
Similar to the BZ mechanism, the magnetic coupling ef- 
fect also extracts rotational energy from the spinning BH. 
Only if the BH spin is initially small would the magnetic 
coupling act as an additional spin-up process. A similar 
discussion on this aspect was made by Dai & Liu (2012) 
within the context of the magnetar central engine model. 
In more general cases, the magnetic coupling effect would 
not significantly affect the BH spin evolution (Lei et al. 
2009). 
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